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We study the finite-temperature optical conductivity a(cu,T) of a granular metal using a simple 
model consisting of a array of spherical metallic grains. It is necessary to include quantum tunneling 
and Coulomb blockade effects to obtain the correct temperature dependence of a, and to consider 
polarization oscillations to obtain the correct frequency dependence. We have therefore generalized 
the Ambegaokar-Eckern-Schon (AES) model for granular metals to obtain an effective field theory 
incorporating the polarization fluctuations of the individual metallic grains. In the absence of 
intergrain tunneling, the classical optical conductivity is determined by polarization oscillations of 
the electrons in the grains, ct(uj) = — (ine 2 fuj/m) / (lo 2 — uj 2 — ijtj|/rg ra i n ), where uj r = e\J (47r/3m)n 
is the resonance frequency, T~ rain is the relaxation rate for electron motion within the grain, and / 
is the volume fraction occupied by the grains. At finite intergrain tunneling, we find that a{to) = 
— (ine 2 ujf/m)/(u> 2 —lo 2 — i|aj|/r ro i) + ctaes(^, T), where is the total relaxation rate that includes 
the intragrain relaxation rate T~^ in as well as intergrain tunneling effects, and <taes(<^',7 1 ) is the 
conductivity of the granular system from the AES model obtained by ignoring polarization modes. 
We calculate the temperature and frequency dependence of the intergrain relaxation time, T(lu, T) = 
t^j 1 — T~* in , and find it is different from <7aes(<^, T). For small values of dimensionless intergrain 

tunneling conductance, g <C 1, the DC conductivity obeys an Arrhenius law, <7aes(0,T) ~ ge~ E °' T , 
whereas the polarization relaxation may even decrease algebraically, T(u>, T) ~ (g / E 2 )[T 2 + (uj /2ty) 2 ], 
when u,T < E c . 



I. INTRODUCTION 

An inhomogeneous mixture of metallic and insulat- 
ing phases exhibits a transition between bulk metallic 
and bulk insulating behavior. When the volume frac- 
tion of metal is large, the composite material is a "dirty 
metal" containing isolated impurities; when the volume 
fraction of metal is very small, it is a "dirty insulator". 
Between these two extremes, there is a third state con- 
sisting of large (~ 100A) metallic regions separated by 
insulating walls. Such systems are called granular met- 
als. Granularity can arise automatically; for instance, 
electronic phase segregation has been directly observed 
in the pseudogap phase of cuprate superconductors^ 
and in two-dimensional electron gases in semiconduc- 
tor heterostructures>£ Granular metals can also be de- 
liberately created by sputtering a metal onto an insulat- 
ing substrate by lithographic deposition of quantum 
dots, or by self-assembly of metal nanoparticles coated 
with organic molecules^ Some of these methods allow 
control of disorder. 

Granular metals are very interesting as their transport 
properties — in particular, the DC conductivity — can- 
not be explained by simple extrapolation from the neigh- 
boring metallic or insulating phases. Another probe of 
the metal-insulator transition is the optical (AC) conduc- 
tivity. In this paper we study the frequency and temper- 
ature dependence of the optical conductivity of granular 
metals. We begin by forming comparative and contex- 



tual links with existing literature on transport in dirty 
metals, dirty insulators, and granular metals themselves. 



A. Dirty metals 

A "dirty metal" consists of impurities embedded in a 
metallic host. The electronic states at the Fermi energy 
are delocalized throughout the solid, giving a finite con- 
ductivity at zero temperature. Thermal excitations are 
detrimental to charge transport, so the DC conductivity 
has a "metallic" temperature dependence (da/dT < 0). 
At very low temperatures, electron-electron interactions 
and quantum coherencei&^iifiiAi can conspire to give "in- 
sulating" corrections to conductivity (d(Aa)/dT > 0), 
which are usually weak. 

The optical conductivity is well described by Drude 
theory, 

2 

O-DrudcM = — — 1 , (1) 

m 1 + l\Ul\TDrude 

where n is the conduction electron density and TQrude 
is the relaxation time, which may be temperature- 
dependent. At high frequencies the optical conductivity 
is dominated by electronic inertia, (TDrudc(^) ~ ne 2 /irnui. 
There are small coherence corrections to the Drude result 
at low temperatures. 



2 



B. Dirty insulators 

A "dirty insulator" or "dirty semiconductor" consists 
of impurities embedded in an insulating host. There is 
a finite density of states at the Fermi energy due to im- 
purity states, but these states are all localized, so the 
DC conductivity is zero at T — 0. Conduction occurs by 
thermally activated hopping between bound states, so 
the conductivity has an "insulating" temperature depen- 
dence (da/dT > 0); it obeys a variable-range- hopping 
law of the Mott 12 or Efros-Shklovskiii^ kind depend- 
ing on whether the long-range Coulomb interaction is 
screened. In this paper we will not be studying the ef- 
fects of long-range Coulomb interaction, and therefore, 
we discuss below only the Mott case. For the sake of 
completeness, a discussion of the Efros-Shklovskii case is 
provided in Appendix lAl 

Mot1*±i showed that the main contribution to optical 
conductivity comes from resonant absorption by pairs of 
states, one of which is occupied and the other empty. 
Mott's argument, which we recapitulate briefly, is valid 
when electron correlations due to long-range Coulomb 
interactions can be disregarded. Let the two states in 
a pair have energies and e^. The resonance condition 
is satisfied when lu = e,j — q. The transition rate Pjj 
in presence of an electric field Ecos(cut) is given by the 
Fermi Golden Rule, 

Pij = 7re 2 V\x ij \ 2 E 2 p(e j ), (2) 

where p(e) is the density of (impurity band) states per 
unit volume and V is the volume of the system. The con- 
ductivity o(lu) is then found by multiplying by lu/^E 2 , 
averaging over all occupied initial states with energies in 
the interval — a-> and e f , and averaging over all unoc- 
cupied final states j. The result isi^ 

Re(aH) = 2ne 2 VLo f de \x ij \l v p{e)p{e + lu). (3) 

The best scenario for a hopping transition between the 
two localized states at low frequencies is that they are 
degenerate and the splitting of the levels to to tunneling, 
Aetj ~ we~ Xij ^ loc is smaller than lu, or in other words, 
the distance Xij between the localized states should be 
large enough: x^ > r u — £/ oc hx(w/io). Here w is an en- 
ergy scale of the order of the relaxation rate^ Localized 
states in a "shell" of thickness £; oc around r u will also 
satisfy the condition for resonance. Using this in Eq.©, 
we arrive at Mott's optical conductivity for a disordered 
insulator, 

Re(aH) ~ 2ne 2 n 2 m ^/5) 2 (ri-^ loc )rl 

« 2ne 2 n 2 mp (u;/Snt + c 2 ln d+1 (w/cu), (4) 

where 7ii mp is the number density of localized states, d 
is the dimensionality, and 1/5 is the density of states at 
an impurity site. An important assumption in obtaining 



Eq.Q is that there is no inelastic scattering during the 
hopping process. 

At high frequencies, u> 7$> w, the electrons are not lo- 
calized, and the optical conductivity reverts to the Drude 
expression, Eq.Q, with rtj mp as the conduction electron 
density. At some intermediate frequency, the optical con- 
ductivity has a maximum; however, this maximum is just 
due to a crossover between different behaviors, and is not 
associated with any special resonance. 



C. Granular metals 

A granular metal consists of metallic grains embedded 
in an insulating host. The electrons are localized within 
each grain due to the Coulomb blockade. Conduction oc- 
curs by intergrain tunneling of thermally excited charges, 
so the DC conductivity has an insulating temperature de- 
pendence (da/dT > 0). However, a granular metal differs 
from a dirty insulator in that there is a large number of 
states N = (a#/i?)~ d on each grain, so the mean level 
spacing 8 ~ £f/N is very small. For temperatures (or fre- 
quencies) higher than S, these closely-spaced levels may 
be treated as a continuum leading to incoherent or dissi- 
pative transport phenomenaii§iiLi£ii2i22i2ii22i22i2i Inelas- 
tic cotunneling, in particular, is the core of the variable- 
range-cotunneling mechanism of charge transport in a 
disordered granular metalf 2 ^* 2 ^ and has an even greater 
effect on heat transport 2 -^. The low-energy particle-hole 
excitations within each grain also give rise to a metallic 
linear-in-T specific heat. 

The standard model for studying dissipative trans- 
port in granular superconductors was obtained by Am- 
begaokar, Eckern, and Schon (AES) in 1982. 16 This 
model has also been widely used to study normal gran- 
ular metalsi 1 Ti 1 ?i 1 ?i 2 9i 21 i 22 i 2 ?i 24 It describes the competi- 
tion between incoherent intergrain tunneling (character- 
ized by the dimensionless intergrain conductance g) that 
tends to delocalize charge, and Coulomb blockade (char- 
acterized by the charging energy of the grain, E c ) that 
suppresses intergrain tunneling. These are quantum ef- 
fects that are beyond the realm of classical electrody- 
namics and circuit theory. The AES approach is valid 
at temperatures larger than both the mean level spac- 
ing in a grain 8 and the Thouless energy of intergrain 
diffusion^Si 2 ^ In this regime, intergrain transport is in- 
coherent and quantum interference effects are unimpor- 
tant. 

We now turn to optical conductivity. The study of 
optical properties of metal particles has a long history 
and occupies a large body of literaturei^^^^SiSLSSi^ 
Effective-medium theories are perhaps the most common 
approaches^ The earliest of these is due to Maxwell Gar- 
nett who in 1904 proposed using frequency-dependent 
dielectric functions in the expression for the effective di- 
electric constant of the granular metal that had been ob- 
tained from electrostatics^ Thus if e m (io) and £i(w) are 
the bulk dielectric functions of the metallic and insulat- 
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ing phases, and e e ff(w) is the effective dielectric constant 
of the composite, 



£ c ffM - £»M _ j. e m {u) - £i(w) 



e e ff(w) + 2£i(w) 



l (w)+2e i (w)' 



(5) 



where / is the volume fraction of the metal. Alterna- 
tively, following Bruggeman^& one can treat the granu- 
lar system as a fraction / of metal and 1 — / of insulator 
immersed in an effective medium. The effective dielectric 
function is obtained by solving 



f- 



n {u>) - £ e ff(u) 

+ 2e cff (cj) 



_ Sjjij) -e eff H = 
Ei(w) + 2e cff H 



(G) 



In 1908 Mie recognized the importance of polarization 
oscillations for the optical conductivity^ The classical 
optical conductivity of a clean spherical metallic grain 
can be inferred from the equation of motion of the elec- 
trons. Suppose an external field E ext e luJt acts on a spher- 
ical metallic particle and induces a polarization P. From 
classical electrodynamics, the field ~Ei n t inside the par- 
ticle is ~Emt — E ea;t — (47r/3)P. Using the equation of 
motion of the electrons, — w 2 x w = — eEj n j/m, together 
with the definition of the current density j w = —ineojyL u 
and its relation to the polarization, j w = iojP, and the 
external electric field, j w (q = 0) = a(oj)E ext , we arrive 
at 



a(u>) 



ine 2 f 



m 
At: ne 2 
3 m 



(7) 
(8) 



u> r , the frequency of resonant polarization oscillations, is 
smaller than the plasma frequency of the bulk metal, 
ey/ Aim/m, by a factor of 1/V3, and depends on the 
shape of the grain but not on its size4^ At very high 
frequencies, u> 3> oj r , the optical conductivity approaches 
that of a free particle, because the inertia of the electrons 
prevents them from screening the external electric field. 
If the electrons in the grain have a finite relaxation time, 
the equation of motion (— uj 2 + i|w|/r re i)x w = —eEmt/m 
gives 



me 



CO 



m uj 2 — ui 2 — i\u>\/T Te \ 



(9) 



The Mie approach is entirely classical. In order to 
capture the temperature dependence of a, which is deter- 
mined by tunneling and charging effects, one has to use 
a quantum treatment such as AES effective-field theory. 
In the original AES model, the Coulomb interaction is 
approximated by a capacitance matrix; the electrostatic 
potential is uniform on each grain (although it may fluc- 
tuate in time). This amounts to assuming that the elec- 
trons are massless and can instantaneously redistribute 
to suppress potential variations within the grain. Such 
a "monopole" approximation is adequate insofar as DC 



(a) 




(b) 







AE~e2 x 2 /R3 



Figure 1: (a) Polarization in the standard AES model in- 
volves charge asymmetry between different grains. For weak 
intergrain tunneling, the energy AE associated with the po- 
larization is of the order of the charging energy of the grain, 
e 2 /R. (b) Polarization due to uneven distribution of charge 
within a grain. Such polarization excitations cost less energy 
than case (a) but are not considered in the standard AES 
approach. The filled large circles denote grains, and e and h 
denote electron excess and deficit, respectively. 



transport is concerned, because the bottleneck in trans- 
port is intergrain tunneling rather than electronic inertia. 
Optical properties, however, depend crucially upon the 
finite mass of the electrons and the possible polarization 
of individual grains (see FigQ. Indeed, a calculation of 
<j(oj) from the AES action alone misses the polarization 
resonance peak completely, and thus severely violates the 
sum rule. 



D. Purpose of this paper and results 

In this paper, we generalize the Ambegaokar-Eckern- 
Schon (AES) model for a regular array of spherical grains 
to include dipole (polarization) as well as monopole 
(charge) degrees of freedom. Using this new effective 
field theory, we are able to calculate the conductivity as 
a function of temperature as well as of frequency: 

1. Using a Kubo formula, we find that the optical con- 
ductivity of isolated grains is mainly due to intra- 
grain dipole oscillations, 



a(w,T) 



ine 2 f 



i\u>\/T ( 



grain 



where T~ ain is the relaxation rate for intragrain 
scattering and consists of, apart from the classical 
Drude relaxation in a bulk metal, additional finite- 
volume effects such as Landau dampingi^iiSiiS 

2. At finite intergrain tunneling, we find that there 
is a small additional "monopole" contribution 



caes(w,T) due to intergrain charge oscillations, 
and that intergrain tunneling also imparts an extra 
width r to the dipole resonance: 



ine 2 f 



i\u\/T„ 



where T Te l 



• grain +r(o;, T). At finite temperature, 
caes(0,T) is finite and gives the DC conductivity 
of the granular array. T(u>, T) depends on the inter- 
grain dimensionless conductance, g, and the grain 
charging E c . It is independent of T~.* in , and has a 
different temperature dependence. 

3. The temperature and frequency dependence of 
the resonance width T(oj,T) is different from 
caes(w,T), especially when u>,T are smaller than 
the effective charging energy of the grains. At large 
u>, T, both r and ctaes become independent of u>, T 
and are proportional to the dimensionless inter- 
grain tunneling conductance, g. Figl^illustrates the 
physical difference between the two. The qualita- 
tive difference in the manner in which intergrain 
tunneling affects ctaes and T cannot be explained 
by a simple effective medium approximation^ 

4. The optical conductivity of a granular metal is 
physically different from a dirty insulator with 
Coulomb interaction even though both show sim- 
ilar features. For both systems, cr(u>) vanishes as 
a power law at low and high frequencies and has 
a maximum at intermediate frequencies. However, 
for a granular metal this maximum is due to res- 
onant polarization oscillations, whereas for a dirty 
insulator the maximum is just due to a crossover 
and is not associated with any resonance. 

Our theory neglects the interactions between dipole de- 
grees of freedom on different grains. However, at high fre- 
quencies (so that the metal dielectric function approaches 
unity) or for small /, our approximation approaches the 
Maxwell Garnett result, Eq.JSJ), when we use the stan- 
dard relation Hea(cj) = —^-lm.e e s(w). In many exper- 
imental situations, such as two-dimensional granular ar- 
rays, it is possible to screen out long-range Coulomb in- 
teraction within the sample by using a gate electrode, in 
which case an effective medium treatment is not neces- 
sary. In any case, the dipole-dipole interactions can be 
modeled with a matrix if necessary, just as the monopolc- 
monopole interactions are included in the original AES 
model as the capacitance matrix. For simplicity, we also 
ignore possible effects arising from the non-uniformity of 
the shape of the metallic particles. 



E. Arrangement of the paper 

In Sec^we introduce our model of granular metals: 
an array of spherical metallic grains with interacting elec- 
trons and a finite intergrain hopping. We then make a 



(a) 





(b) 



Figure 2: Physically different mechanisms involving inter- 
grain tunneling processes that determine the DC conductiv- 
ity a (that measures the escape rate of an electron from the 
grain) and the relaxation rate of polarization oscillations Y . 
(a) Coulomb blockade of intergrain tunneling, especially at 
small values of g, suppresses electron escape from the grains 
leading to a ~ exp(— E c /T) at low temperatures, (b) A polar- 
ization oscillation does not result in a net transfer of charge 
and therefore no strong Coulomb blockade at low temper- 
atures. High temperatures and strong intergrain tunneling 
both wash out Coulomb blockade effects, whereby the two 
processes show the same temperature dependence. Drude the- 
ory predicts the same temperature dependence for a and T. 



multipole expansion of the potential on the grain in terms 
of spherical harmonics up to the I — 1 (dipole) compo- 
nent. In Sec lllII we develop an effective field theory in 
terms of the monopole and dipole components of the po- 
tential fluctuations: this is a generalization of the AES 
theory of transport in granular metals. Optical conduc- 
tivity is calculated in Scc llVI First we consider isolated 
grains and calculate the optical conductivity using the 
Kubo formula approach as well as from the dielectric 
function. The calculation with the dielectric function is 
much less tedious. The result agrees with classical ex- 
pressions for the optical conductivity of isolated grains. 
Next we consider the case of finite intergrain tunneling 
and study the differences from the classical optical con- 
ductivity. An explicit expression for the broadening of 
the polarization resonance, T, is obtained. Its temper- 
ature and frequency dependence are found to be differ- 
ent from the conductivity of the granular metal obtained 
from the AES model. The paper concludes (Sec[Vj with 
a discussion of the results and open problems for further 
study. 
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II. MODEL 

Wc consider the following action for the granular metal 
array, 

+ E/ ^rxJ^+e(^ V x i )]^rx i + 

j «/ TXj 

+ E/ '^•-'•-/•-r (10) 
(ij> Jtx ' x J 

where p(xi,r) = (^jl x .%p TX . T — Q i)O{\xi\—R) is the excess 
electronic charge density at position Xi in the i th grain 
of radius R, £(— iVi) = e(— «V Xi ) — /x = p^./2m — /i, a 
is a lattice translation vector, and £ X1 , XJ is the intergrain 
hopping amplitude. Integrals over r are understood to go 
from to p. We assume the intergrain hopping amplitude 
has a white-noise distribution, 

(WjWi) = l < | 2(5 ( x i - x i)^( x j - x k)- (11) 

where angle brackets denote disorder averaging. 

Next we decouple the Coulomb interaction in Ea. lfTU|) 
through a Hubbard-Stratonovich field, (xj), that has 
the physical meaning of the electrostatic potential. The 
interaction part of the action is 



-^lE / C(x i ,x j )T/ i (x i ,r)y j (x j ,r) 

Ze ij • / TX,Xj 
* J TX; 



(12) 



where J x C(xi, xj) | x .J L Xk = <5(xi — Xk) subject to appro- 
priate boundary conditions at the metallic grains; thus 
C(xi,xj) is proportional to the Laplace operator. 

For simplicity, we will consider grains sufficiently far 
apart so that the mutual interaction of electrons on dif- 
ferent grains is small compared to the interaction of elec- 
trons within individual grains. With this simplification, 
the interaction part of the action becomes 



8 ^E/ Tx Ei( x ,-)-E i( x,.) + 
+ E / ^UVi^r)^, (13) 

j J rxi 

Ei(x, r) = — VVi(x, r) is the electric field at x due to 
charge on an isolated grain at i. The potential away from 
the boundary may be expanded in a basis of eigenfunc- 
tions of the Laplace equation, 



I ill 



V,(x,r) = E 5 i iT » 



R 



lm\ 
1+1 



Yu, 



<R), 



(r > R), 



where r = |x|. Continuity of V[(x, t) at the boundary 
requires that A l ™(r) = Bj (r). For the purposes of this 
paper, it is sufficient to retain just the monopole compo- 
nent (average potential) 

Vi(l = 0; x,r) = | Fio(r)( ^° / ( r T ) ) ; r> R, 
and the dipole components (electric field), 



^ a) (l = l;x,r) = 



vS 1 a) ( T )(R 2 x a /r 3 ), r>R. 



Vn\r)(x a /R), r<R 



Using the definition Ei(x, r) = — VVj(x, r) the interac- 
tion part of the action, Ea. (|13|l takes the form 



R 



2e 2 ^ 

i 

E / *U 



(^o(r)) 2 +E^ii a) W) 2 



V i0 + J2Vn\^/R) 



T/W (W) 



Now make the gauge transformations 



V, TXi Vvxi e - i *"o(r)-iE a y, ( 1 0!, (r)W/ii) ) ( 15 ) 

y i0 (r) = id M V^\t) = idripfiK (16) 

to eliminate V5o(t) and replace Vh(t) by a time- 
dependent vector potential, 



: J TX; 



iO" 



VVx, 



E/ ^x 1 [9 T +^HV Xl - ¥ ) i l/i?)]Vrx 1 . (17) 



The gauge transformations also dress the tunneling cle- 
ment in Eg. Ijl L)| l with monopole (I — 0) and dipole (1 = 1) 
phase fluctuations, 



t, 



x exp 



tx 1 ,x J (r)=t xl)Xj .e i ^W-^W)x 



(i/R)J2(<P^W-<p^\r)x? 



(18) 



III. EFFECTIVE FIELD THEORY 

Integrating out the conduction electrons results in an 
effective action for the I = and / = 1 phase fluctuations: 

Ses[V>] = 2^2 E/ [(^^ i0 ) 2 + I^Vill 2 

-tr In [G-i <5ij - (i xi , xj $ J)1+B + i «-» i + a)] , (19) 
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where 



^=^ + ^(p Xi -^j -A. (20) 



R 



is the inverse of the electron Green function on grain i 
in the absence of intergrain tunneling, and t X; Xj is the 
dressed tunneling amplitude defined in Eo. ((T%)) and the 
bare tunneling i Xi .Xj has a gaussian distribution as in 
Ea. (|ll|l . We study first the effective field theory for iso- 
lated grains and then consider the effect of finite inter- 
grain tunneling. 



A. Isolated grains 

In the absence of tunneling, the "bare" effective ac- 
tion [</?] is obtained by expanding the determinant in 
Eq.^HJ up to second order in <£ii(r), 



2e 2 



E J [(<^io) 2 + \d T <fii\ 2 



+ -/ V,i(r)-p 1 , 1 G!° , (x,,xi ; T,T') 

m J rT'xjXi' 



(21) 



where g[°^ — — [d T + ^Px ; — A 4 ] * i s tne bare electron 
Green function, 

G (o) (x .^ ;Ti / ) = rV^e-"("') 



A , n 



^GS ) (r,r')^(x!)V^(x i ). (22) 



Note that 



-iu n (t— t') 



A n ^"^ " ?1 

= -2ii/(e F )r^e- w " (T - r,) cot" 1 j/„t c 

(23) 

where r c ~ e^ 1 is a short time cutoff. For |t — t'\ 3> r c , 
this simplifies to 



(0)/ is _ irTv(e F ) 



sin7rT(r — r') 



(24) 



We shall use this expression unless stated otherwise. 
Ea. H2f|) can be presented in a more recognizable form 



■ a J t.t' 



4tt i? 
T2? 



(25) 



in terms of the bare electromagnetic response function of 
the \ th grain, 

Kf(r - t 1 ) = 8 a P^-J^ Gf ( Xi) x,; r, t>(t _ /)+ 

(26) 

where V = (47r/3)i? 3 is the volume of the grain. In 
a bulk metal, the two terms in the electromagnetic re- 
sponse function in Eq. i|26[l would correspond to the dia- 
magnetic and paramagnetic parts of the bulk conductiv- 
ity <j(oj) = K(u)/ioj. In a finite system, the situation is 
trickier. The frequency dependence of K is 



Kf{iu> m ) = 5 at) 



ne 
m 



2 /(6a)6,aa'1p£av1 2 



mN 



XX' 



XX' 



(27) 



where A, A' label the eigenvalues of the free electron 
Hamiltonian of a grain, £,aa' = £,ix - fiA', /(£) = 
[e^ + is the Fermi-Dirac distribution function, and 
pfxx' — (Mp?\^')- 71 — N/V is the number density of 
electrons in a grain. If the temperature (or frequency) is 
much smaller than the level separation S = ~ cf/N, 
we expand the right hand side of Ea. l|27() in ascending 
powers of ui m : 



ne 
rn 

,2 \ ^ 



2 v /(6a)K aa ,| 2 

mN ^i.AA' 



f(ZixMx\>\ 



mN 



AA' Sl > AA ' 



+ 0(<4,). (28) 



The static part of Eq. I|28[l can be shown to vanish using 
the Reiche-Thomas-Kuhn sum rule4^^ 



Kaa-I 

6, 



= -l, 



AA' 



(29) 



along with the identity f{£,ix) — N. Combining 
Ea. (|25|) and Eq. l|28l) . one finds that the surviving contri- 
bution in Eq. (|28() makes a finite size quantum correction 
to the RPA dielectric constant 28 i 29 . 30 i 31 . 32 i 33 , e RPA = 

2e 2 v-v S^)\Pi.xy 

2^aa' ?3 ; 



I 



1 — (kFao){R/ao) 2 , where 



ao is a small length of the order of a lattice constant. 
As a result, even for metallic grains a few tens of lattice 
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constants across, the static dielectric constant rapidly ap- 
proaches bulk values (where is it infinity) , and the polar- 
izability, a = R 3 (e RPA - l)/(e RPA + 2), approaches the 
classical value, a c i assica i = R 3 . 

The sum rule enables us to recast the electromagnetic 
response function as 



2e 2 2V mxM^l 



(30) 



which is a known result. For the rest of the paper, unless 
stated otherwise, we shall assume that the temperature 
(or frequency) is much larger than the level separation 
(T/S > 1). Then, using Eq.® and Eq.JEJ, we obtain 



ra/3 



ne 



T/S > 1, 



(31) 



that is, the value for the clean bulk metal; this is the 
diamagnetic response due to electron acceleration in an 
electric field. Hence, Eci. (|25|l becomes 



(32) 

where u r is the resonance frequency for a metallic sphere, 
9 47r ne 2 



3 TO 



(33) 



we introduced in Eq.©. In a collisionlcss bulk metal, the 
paramagnetic part of the electromagnetic response func- 
tion defined in Ea. l|26|l vanishes. However, in a finite- 
size grain, if one can treat the quasiparticle excitations 
in the grain as a continuum (this is so if the tempera- 
ture is not too low, T 3> i5), the paramagnetic part is 
finite and gives rise to a finite relaxation of the oscilla- 
tions through disintegration into incoherent particle-hole 
excitations^i*^ The relaxation time has been shown in 
numerous works^&iSi to be of the order of the time of 
flight, R/vf- Physically, the relaxation is due to Landau 
damping of plasma oscillations at a finite wavevector: the 
minimum wavevector in a grain of size R is of the order of 
tt/R. Other inelastic processes such as phonon scattering 
will also contribute to relaxation. 



B. Finite intergrain tunneling 

We now obtain the effective field theory when 

intergrain tunneling is finite. At not too low 

temperaturesj22i2i T ^> max (|i| 2 <5, 5), and for large 

enough 43 grains (kpR) 2 1, it suffices to expand the 



electron determinant in Eo. i|19|) up to 0(t 2 ), 
S^>] = / ^.x; +a (T% i+a , Xi (T) 

i,a ^T'xixlxi+.x:^ 

x G Vil (xi, x ; ; r, r )G Vl+tll (x i+a , Xi+ a ; r , t) 

= E / G(x ; , x ; ; r, t')G Vi+a1 (x i+a , x i+a ; r', t) 

: ~ rr ' x i x i + a 



s *(mo( T ')-¥>ij,o( T )) 

Jil/RKipuiT^-ipnir))-^ „-i(l/il)(Vji(T')-V»ji(T))- X j 



(34) 



where <Pij,o( T ) = ¥'io( T ) ~ Vjo (t). The cpn dependence 
in Eg . I|34[l comes from the exponential as well as from 
the the Green functions, G Vil (xi, Xi; r, r'), etc. We show 
in Appendix [5] that the contribution arising from the 
expansion of G Vil (xj, Xi; r, r') etc. in powers of ipn is 
insignificant compared to that coming from the expo- 
nential. Therefore we expand only the exponential in 
Ea . (|34|l up to second order in tpn etc. and ignore the tpn 
dependence of the Green functions. Thus the tunneling 
part of the effective action is 

S*™^| 2 E / ni,i + a(r,r')GW(r ) rOGW iA ,(r' 1 r) 

i,a;AA' 
i,a;AA' rr 

x[r 2 A ( V ii(r')-¥>ii(r)) 2 

+ »f+a,A'(*'i+a,i(r / ) -¥>i+a,i(r)) 2 ] (35) 

where 

ni.i+a(T, T') = COs(( i 9i ! i +a ^o(r)-^i4 +a ,o(T')) (36) 

and rf x = (iA|f 2 |iA) are the matrix elements of f 2 = 
|x| 2 with the eigenstates of grain i. For a spheri- 
cal grain, the eigenfunctions are spherical Bessel func- 
tions jn{n n ir I R)Yi m (8 \4>), where n n i is the I th zero of 
j n (x). Numerically evaluating the matrix elements we 
find they range between 0.28i? 2 and R 2 . In particular, 
lim„^ 00 (r 2 )„i = R 2 , and lim/^ 00 (r 2 )„ ; = R 2 /'i. Thus 
the matrix elements of r 2 do not vary strongly and are 
of the order of R 2 . Hence Eq. (j35(l may be written 



vgT 2 ^ f n i4+a (r,rQ 
2 Jtt> sin 2 itT(t - r') 

i,a v ' 

+ ngT 2 by, r n i , i+a (r,r / ) 
2 <H Jtt> sin 2 itT{t - t') 

i,a v ' 

x [{<Pu{t') - Vh{t)) 2 + (Vi+a,i(r 7 ) - ¥>i+a,i(r)) 2 ] ; 

(37) 

here b ~ 0.1 is a constant, g is the dimensionless inter- 
grain tunneling conductance 



g = 2ir\t\ 2 u(e F ) 2 , 



(38) 
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Eas. l|32[) and l|37|) form the effective action, S'cffM — 
S^g [ip] + Slf[ n [tp] , which generalizes the AES action to 
include the physics of dipolar oscillations. This may be 
presented as S e s[ip] « Saes[</?o] + S P oi[tpo] where 



T 

2(e 2 /R) f 



T 



y^r i;i+a (^ m )|cj m |x 



4(e 2 /i?) 

[<Pil(u m ) ■ tfiil(-LJ m ) +¥>i+a,l(^m) ' ¥>i +a ,l(-W m )], (39) 

where 



5*AESM = 2(e2/fl) Z J^rVio) 



E 



ni,i+a(T,TQ 

rT < sin 2 ttT(t — t') 



(40) 



is the standard Ambegaokar-Eckern-Schon (AES) model 
for normal granular metals and 



sin 2 (ttTt) 



(41) 



The quantities II and T are functionals of </?o- I n princi- 
ple, fluctuations of ipo and of ipi can influence each other 
since they both appear in S po \. In practice, it is suffi- 
cient to calculate the correlator (II[y>o]) f° r <5aes alone, 
and to use this mean value in S po \ to determine the fluc- 
tuations of tpi. To justify this, we show that fluctuations 
of (fx have a negligible effect on the "kernel" for ip . If 
in Ea. H37fl we average over the fields <pu using their bare 
propagator in the absence of tunneling (see Ea.l|3*2jl). 



<(^ii(r)- V ii(0)) 2 )=3( e 2 /i?)T^ 

n 

,e 2 /R 



(1 - e - l ^ T ) 



coth(w r /2T)(l - e-^l T l). 

2ui r 

(42) 



Thus at long times, w r |r| ^> 1, the correction to the 
tunneling term of the AES model due to dipole modes is 
smaller than the bare value by a factor of (e 2 / 'R)/ui r . In 
most common cases of granular metals, this ratio is of the 
order of 1(T 2 , as E c ~ 10 2 K and u r ~ 10 4 K. At short 
times, w r |r| -C 1, the correction is smaller than the bare 
value by a factor w r |r| <C 1. Thus under most common 
physical circumstances, our approximation is valid. 

For finite tunneling, the propagator for the dipole 
modes is that of a damped harmonic oscillator, 



^<5jj(e 2 /i?) 

ijj 2 , + uj 2 + r(uj n )\u n \ ' 



(43) 



where the resonance linewidth is 



+a- 



(44) 



IV. OPTICAL CONDUCTIVITY 

In this section we calculate the optical conductivity of 
isolated metallic grains and then generalize it to finite 
intergrain tunneling. For isolated grains, we show that 
the optical conductivity may be obtained in two ways: 
directly from the Kubo formula, and from the dielectric 
function. 



A. Isolated grains 

1. Kubo formula approach 

We first calculate the optical conductivity for isolated 
grains using the Kubo formula approach. For this we 
introduce an infinitesimal vector potential A rx that cou- 
ples to the current j TX and is related to the electric field 
through E rx = i<9 T A rx . The electronic kinetic energy 
becomes 



e(prxi) -* e(p Xi - f A TXi ) 

= e (pxi) — ^;A TX . • p x . + A TX ., 



(45) 



where p Xi = — iV Xi and p TXi = p Xi — ^(Ph(t), and we 
have chosen the gauge V- A = 0. The optical conductivity 
tensor <r Q/3 is the coefficient relating the q = component 
of the current, 

j?[A;qfl =-±J dxe-^J ^fields) j^e'*™, 

(46) 



to the q = component of the electric field E@ x 
id r A0 



a a f 3 (T,r') = Tc 2 J2^(^ m y 



«#[A;q = 0] 



— ZUJjn (t — t) 



M£(q = 0) 



(47) 



A=0 



Here Zq = Z[A = 0]. Analytically continuing to real 
frequencies gives the well-known Kubo formula for the 
optical conductivity, 



ct q/3 (w,T) 



IC 



ujZqV 
I> (fields) 



5 2 S[A] 



<55 [A] 6S[A] 



SA^, 5A% 



(48) 
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where V is the volume of the system. The first term in 
Eq.|@8J, as we shall see below, represents the inertial re- 
sponse of the electrons in the bulk metal. The second 
term, which vanishes in the bulk, makes a finite con- 
tribution in the granular metal. We denote these two 
contributions as 

<t*V T) = <f rtii >, T) + <f itc R ( W , T). (49) 

We can show that the "inertial" term is 



inertial ( W > T ) 



u-y 



ujmV 

2 



E 



ie 



ujm 2 V 



2i; E 



J cfoq (G vn (x ; , x ; ; r, r)) - 



x G Vil (x-, x ; ; t, 0)G Vil (x;, x£; 0, r)) 

(50) 

This can be expressed in terms of the response func- 
tion we defined in Ea. (|26|) . For simplicity we as- 
sume that all grains in the system are identical. Also, 
as in Appendix [El we approximate the Green functions 
G Vil (x-, xi; r, 0) etc. by their bare values. Then 

,V,2 



aB 

T 

inertial 



if 



xG[° ) (r,0)G I ^(0,r)(^ 1 (r)^ 1 (0)) 



-,(0) 



2 R 2 V 



E 



(51) 

where / is the volume fraction occupied by the metallic 
spheres. 

The second term on the right hand side of 
Ea. (|51fl is smaller than the first by a factor of 
i/(eF)(e 2 1 'R) j ' (Nuj r mR 2 ), where N is the number of con- 
duction electrons in a grain. Since l/(mR 2 ) ~ S ~ ep/N, 
and v{ep) ~ N/cf, the second term is smaller by a factor 
of about l/N. This is a small number since the number 
of conduction electrons in a grain in typical systems is 
of the order of 10 4 . We have, dropping this term from 
Eq.(EU, 

if 



a/3 

7 

inertial 



ifne 2 
mw 



(5 Q/3 , 



(52) 



where we used Ea. (|31|) in the second line. This is indeed 
of the form of an inductive contribution. 

Now consider the finite size contribution to the con- 
ductivity described in Eq.lgS) and Ea. lj39"|) . 

•,2 



aB 
finite 



ie 



torn 



dxjdxj 



x (Prxi-Poxj^ii ( x ii x i! T : T ) x 

x G Wl (xj,Xj;0,0)) 



(53) 



The diagonal matrix elements of the momenta p are iden- 
tically zero in a finite system, (A|p|A) = 0, therefore we 
discard in E!a. (|5Hfl terms of the type 

/dxip^G^ (xi.Xi ;t,t) =0; 
this simplifies the finite size contribution to 



aB 

T r 

finite 



r(^) 



ie 



ujm 2 R 2 V ' 

ij;AA' 



xGg ) (r,r)Gf J ^(0,0)K 1 (r)^ 1 (0)) 



■,(o) 



(54) 



Here we have as usual approximated the Green functions 
G(p n (x^ Xi; r, r) by the bare values. Evaluating Ea . f54p 
gives the following finite size contribution, 



ff WteR( W » T ) = 



ine 2 JV/ (e 2 /i?) 



u)m 2 R 2 to 2 
ine 2 f 



U7. 



uim — (uj + i0 + ) 2 + lo 2 . 



(55) 



Here N = (47r/3)ni? 3 is the total number of conduction 
electrons on a grain, n is the conduction electron density, 
and we used to 2 = (4ir/3)ne 2 /m. Adding the inertial and 
finite size contributions from Eq.JS^J) and Ea.(|55|). we 
arrive at the optical conductivity for isolated spherical 
grains, 



o- a8 (oj,T) 



ine 2 f 



m lu 2 — uj 2 . — i\uj\0 + 



(56) 



Eq. (|56|l agrees with the expression for the optical conduc- 
tivity in Eq. Q that was obtained from a simple analysis 
of the equation of motion of the electrons in a clean grain. 
For a finite intra-grain relaxation time, r gra in, the optical 
conductivity takes the form 



a al3 (uj,T) 



ine 2 f 



IUJ T, 



(57) 



grain 



2. Conductivity from the dielectric function 

The optical conductivity of isolated grains that we ob- 
tained from a tedious Kubo approach could also be in- 
ferred from the dielectric function. In a gaussian theory, 
the dielectric function e for a single grain can be ex- 
tracted from the effective action, 



S = -S~2 E / ^(r-r'-.Eil.r.r,;,/-.. m 



where E£(t) are the multipole components of the elec- 
tric field (see Ea. (|13|i fFl at the grains. The lowest possi- 
ble angular momentum component of an excitation on an 
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isolated grain is I = 1. That is, in the absence of inter- 
grain tunneling, the simplest response to an electric field 
is a uniform polarization. Thus we need to consider only 
E9[{t) = Vft/R = idr^i^/R. Furthermore, because of 
the high energy io r ~ e^J n/m associated with the dipole 
excitations, we can safely neglect in the effective action 
terms with higher powers of <pn(r). 

The following relation can be gathered from Ea. (|32[) . 
Eq.(G2J and Eq.©, 

S al3 ine 2 f 8 a0 Lu 



T a/3 



9=0 



ine f 
mio £ii(w,T) 



- "7M 
(59) 

Such a cross relation has been discussed, for instance, 
by Hopfield^i in 1965. Physically, the imaginary part of 
the dielectric function is associated with relaxation, so a 
stronger relaxation implies weaker conduction. 



The conductivity due to the monopole part has been ob- 
tained elsewhereSi in the context of the AES model, 



af(cv,T) 



AESV J\n ^-iu> 



(61) 



where a is the intergrain distance. It consists of diamag- 
netic and paramagnetic parts, 



K^f a ( T )=S^e 2 ngT 2 / (6(r) - S(r' - r))x 



sin (ttTt 



-(cOS(^i,i +ect ,o(T) - ^i,i+e Q ,o(T'))), 

(62) 



B. Finite intergrain tunneling 



AES \~) 



-^^(X«(r)Xr(0)), (63) 



The Kubo approach is the most reliable way to cal- 
culate the optical conductivity, but, as illustrated in 
Sec lIV Al it is very tedious even for an isolated sphere. 
At finite intergrain tunneling, an even larger number of 
terms involving both intragrain and intergrain currents 
would have to be calculated. We also showed that for 
gaussian models, the dielectric function could be used 
to obtain the conductivity with significantly less effort. 
However, as the discussion below shows, the theory is 
gaussian only in the two extreme cases of isolated grains, 
g = 0, or strongly coupled grains, g ^> 1. So we resort 
to a combination of the Kubo and dielectric function ap- 
proach, using the Kubo approach for multipole modes 
that cannot be considered in a gaussian approximation, 
and retaining the dielectric function approach for modes 
that are effectively gaussian. 

At finite intergrain tunneling, an electric field can 
cause intergrain polarization (opposite charges on ad- 
jacent grains) as well as intergrain polarization. We 
must therefore consider the contribution of the monopole 
modes (fiofr) in the dielectric response function. Tun- 
neling events are accompanied by fluctuations in electro- 
static energy which can be large, of the order of e 2 /R, 
when intergrain tunneling is weak, g -C 1. Therefore 
for weak but finite tunneling, we must consider non- 
gaussian contributions for the monopole modes, Vio(r) = 
id T ipio(T). This is clear from the effective field theory 
at finite tunneling given by Ea. (|39|) and Ea. l|40l) . On 
the other hand, for strong intergrain tunneling, j > 1, 
monopole fluctuations are small because charges can eas- 
ily flow to neutralize potential differences between the 
grains. In this case, we again have an approximately 
gaussian theory for the Lp- l0 {r) modes. We write the to- 
tal conductivity as a sum of the monopole and dipole 
contributions, 



t q/3 



a/3 



a/3 



( W ,T). 



(60) 



where 



X?(t) = eirgT 2 J 



sin 2 [7rT(r - r') 



x (sin(^ M+eT ,o(r) - ^,i+e T ,o(r / ))>- ( 64 ) 

In order to remind us of the AES origin of the / = 
component of the conductivity, we rename <7o(oj,T) to 
ctaes (w,T). 

The contribution to the conductivity from the dipole 
part is written in terms of the / = 1 component of the 
dielectric function (see discussion above and in the pre- 
vious section), 



af(co,T) = - 



ine 2 f 8 aP 



mcj en(uj,T) 
ine 2 f 8 al3 uj 
m lu 2 — to 2 — iT(u>, T)\u>\ 



(65) 



Here T(lj,T) is the intergrain relaxation rate defined in 
Eq. (l41|) . r involves the cosine correlator n defined in 
Eq. (I36|) , and thus closely resembles the diamagnetic part 
of the AES conductivity, Ea.fity. However, the full AES 
conductivity behaves very differently because of the para- 
magnetic contribution, Eq. H63|) . 

In the presence of intragrain relaxation mechanisms 
such as impurity scattering or boundary scattering, we 
expect that 



ine 2 f S al3 Lj 
m lj 2 — U) 2 — i\uj\/T rc \ 
(66) 



for consistency with Eq.® and Eq. (|57f) . where the to- 
tal relaxation rate t~1 is given by the Matthiessen rule, 



rcl grain 1 
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C. Some special cases 



logarithmically^! with temperature and frequency, 



The final expression for the optical conductivity con- 
tains the AES conductivity oaes^, T) and the resonance 
width due to intergrain tunneling T(u,T). oaes^i T 1 ) = 
<Jq(lu,T) has been explicitly defined in Ea. l|tjl|) through 
Eq.(jni}, and T(uj,T) has been defined in Eq.J2B and 
Eq. (|44H . Below we discuss a few special cases for a regular 
three dimensional array. Throughout we assume that the 
frequency lies in the range lj/T ^> 1, u>t c ~ oj/ef <C 1, 
and the temperature much smaller than the charging en- 
ergy, E c /T » 1. 

(a) Consider first small intergrain tunneling conduc- 
tance, g -C 1. As the calculations are very complicated, 
we refer the reader to Appendices |0 an d El f° r the de- 
tails. We show there that at frequencies much larger than 
the charging energy, the conductivity tends to saturate, 
caes ~ ff(e 2 /a). The same goes for the polarization res- 
onance width, r w Abzg(e 2 / R), where z is the grain co- 
ordination number and we used Ea . (|41fl and Ea . Q44JI . If 
the frequency is much smaller than the charging energy, 
the conductivity is dominated by thermal excitation of 
quasiparticles and obeys an Arrhenius law, 



a 



(67) 



In contrast, the resonance width does not obey an Ar- 
rhenius law: 



Abzg 



e 2 47r 
R 3E 2 



T 2 + (w/2tt) 2 ], uj<^E c 



(68) 



Suppose the charging energy is small compared to the res- 
onance frequency, E c <C to r <C ^f- As significant changes 
in oaes and T occur on the scale of the Coulomb block- 
ade energy E c , the frequencies in the vicinity of the res- 
onance are too large for Coulomb blockade physics to be 
significant. In this case, T sw Abzg 2 (e 2 / R) is practically 
independent of frequency and temperature. 

Consider now the case where charging energy is large 
or comparable with respect to the resonance frequency, 
uj r < E c <C £f- This can happen if the metal has a low 
enough conduction electron density, a large effective mass 
for the electrons, and/or small grains. Increasing the 
volume fraction of the metal is another way in which 
the resonance frequency may be reduced; we shall see 
in Sec II V Dl that to r renormalizes to u>* = u> r y/l — / as 
/ is increased. This regime is very interesting because 
the resonance is in the low frequency regime {ui <C E c ) 
for Coulomb-blockade physics. So near the resonance 
to = to*, while (taes still obeys an Arrhenius law, Ea. H67|l . 
the temperature dependence of T can be qualitatively 
different from caes- One expects here, following Ea. (|68|l . 

t(lu;,t)^(t 2 + (^) 2 )/e 2 . 

(b) Finally consider large intergrain conductance, g 3> 
1. In this case both oaes(w,T) and T{u>, T) evolve 



caes(^,T) ~ g(e /a) 

„2 



IV, T) w Abzg 



R 



1 

ngz 
1 

irgz 



In 



In 



gE c 



max(w, T) 

gE c 

max(w, T) 



down to exponentially low temperatures and frequencies 
when perturbation theory is no longer valid. Below such 
low temperatures, the physics is similar to the g <C 1 case 
discussed above. 



D. Comparison with Drude theory 

According to Drude theory, a bulk metal will have a 
frequency-dependent dielectric function 



EDrudcV) — 1 — 



Anne 



1 



(69) 



which when substituted in the Maxwell-Garnett formula, 
Ec[-©7 yields the effective dielectric function for a homo- 
geneous system of metallic grains in vacuum, 



EeflV) 



to 2 - (47me 2 /3m)(l + 2/) - ioj/Stdt 



ide 



to 2 - (47me 2 /3m)(l - /) - iw/3r Dr ude 



(70) 



From Ea. (|70|) one then infers the optical conductivity 
(Tmg(^) for the granular system in the Maxwell-Garnett 
approximation, 



Re ctmg{u) 



fne 2 



w 2 /(3r Drude ) 



m \u>* 



(w/3TDrudeJ 



(71) 



Ea. (|71|) is not strictly correct because TD ru de does not 
include the Landau damping2i^2i2S that exists in the 
metallic grain but is absent in the bulk. (See also the dis- 
cussion in Sec lIIin Besides, matching Ea. (|7T)> to the cor- 
rect DC conductivity requires that TDrude ~ vaes(T) ~ 
exp(— E c /T), whereas we have shown that the resonance 
width r has a different temperature dependence from 
caes- Evidently, classical arguments are unable to ex- 
plain the full behavior of cr(u>, T). 

Pending a proper theory of long-range interaction of 
dipoles in the granular metal, we nevertheless propose 
that the optical conductivity of the granular metal in the 
Maxwell-Garnett approximation is given by Ea. (j71j) with 
^rain + r replacing (3r D rudc 



' rel 



Re a(to,T) = oaes(w,T)- 
fne 2 



w 2 /r rcl 



m [^-^ r (l-f)] 2 + (u/T Iel y 



(72) 



The Maxwell-Garnett result, Ea. (|72l . agrees with the 
dipole contribution in our stronger result for the opti- 
cal conductivity, Ea. l|66|) . that was derived for a dilute 
granular array, /< 1. Eq.(72J shows that the resonance 
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frequency undergoes an infrared shift, lu* — oj r ^\ — /, 
as the volume fraction of the metal is increased. This 
dependence has been previously obtained^ in the liter- 
ature. One must take care not to extend Eci. (|72|) all the 
way to / = 1 because the validity of our effective field 
theory is limited to the insulating phase of the granular 
metal. 



V. CONCLUSIONS AND DISCUSSION 

We have developed an effective field theory of granu- 
lar metals which is a generalization of the Ambegaokar- 
Eckern-Schon (AES) action to include polarization de- 
grees of freedom. This approach synthesizes the classical 
electrodynamic theories of Maxwell Garnett and Mie and 
the quantum mechanical AES model for dissipative trans- 
port in order to capture both finite- frequency and finite- 
temperature effects. It is valid at temperatures larger 
than the mean level spacing 8 in a grain. 

Using this effective field theory we have calculated 
the frequency and temperature dependence of the opti- 
cal conductivity of an array of spherical metallic grains. 
We have shown that the temperature dependence of the 
polarization resonance width T differs qualitatively from 
that of the DC conductivity for frequencies and tem- 
peratures much smaller than the charging energy of the 
grains. While the DC conductivity obeys an Arrhenius 
law at low temperatures, T decreases only algebraically as 
a function of frequency and temperature. We believe this 
prediction can be tested in experimental situations where 
the condition u) r */l — f < E c can be satisfied. This can 
occur in systems where the conduction electron density 
is low, the effective mass is large, and/or the grains are 
small, and the volume fraction of the metal is large (while 
still remaining in the insulating phase). This qualitative 
difference between the temperature dependences of the 
DC conductivity a(0, T), and the collective mode damp- 
ing r(0,T), obeyed in certain granular metals is quite 
unlike the behavior seeni& in pinned sliding density wave 
compounds where the temperature dependence of the col- 
lective mode damping is the same as the DC conductivity. 
Such a difference could perhaps be used to distinguish 
between granularity arising from spontaneous electronic 
phase segregation in strongly correlated electron systems 
and density wave order. 

To keep our analysis simple, we have, in our field the- 
oretical treatment, ignored electrostatic interactions be- 
tween monopoles (charges) and dipoles (polarizations) on 
different grains. Strictly speaking, this is correct only in 
a dilute granular array (/ <C 1) or at frequencies higher 
than the polarization resonance. Renormalization of the 
resonance frequency due to the presence of neighboring 
grains, even in the absence of tunneling, is one effect that 
is lost in this approximation. Pending a general field the- 
oretical treatment of long-range interaction of dipoles, we 
have used our result for the optical response of a dilute 
array of grains as an input in a Maxwell-Garnett effective 



medium approximation to obtain the optical conductiv- 
ity at larger values of /. The shift that we obtain in the 
resonance frequency as a function of / agrees with earlier 
results in the literature^ 

Another aspect we have not considered is disorder, 
both in intergrain tunneling conductance and as a ran- 
dom background potential due to quenched impurities in 
the insulating part. In presence of strong disorder, the 
DC conductivity obeys a soft-activation law <r(0, T) ~ 

instead of an Arrhenius law; it should be in- 
teresting to consider the effect on optical conductivity. In 
principle it is possible to study the effect of both kinds 
of disorder in our scheme. 

Finally there are some fundamental limitations on 
AES-inspired treatments. Like the AES model, our dis- 
sipative transport model is limited to the insulating side 
of a metal-insulator transition, and cannot describe the 
optical conductivity through the transition; it also ne- 
glects quantum coherence effects, which are important 
at T < 6. 

The present level of rigor in our calculation is insuf- 
ficient to study the various /—sum rules obeyed by the 
optical conductivity^ Our model is justified only for fre- 
quencies much smaller than the bandwidth, u> <C t^ 1 ~ 
tp. At higher frequencies, or in other words for times 
shorter than the cutoff, r <C t c (see discussion following 
Ea.l|2*3)ll. the dissipation kernel, T 2 / sin 2 (7rTr), that ap- 
pears in the tunneling terms in the effective field theory, 
Eci. (|37|) . is no longer valid. 
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Appendix A: EFFECT OF LONG-RANGE 
COULOMB INTERACTION ON OPTICAL 
CONDUCTIVITY OF DIRTY INSULATORS 

Shklovskii and Efros generalized Mott's treat- 
ment to include the effect of long-range Coulomb 
interactions ,isisl The difference is particularly significant 
when the particle-hole Coulomb energy at hopping dis- 
tance r w exceeds the optical frequency: e 2 /sr u 3> u>. Here 
e is the dielectric constant of the medium. Physically, in 
the presence of Coulomb interactions, transitions to the 
final state tj can be made from an occupied level with 

2 

energy in the range ej? — u> ^— < < £f- Modifying 
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the limits in Eq.Q accordingly, 
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Re(aH) ~ 2^e 2 n 2 mp ^ + j- ) (rtT^ior.yr 



(Al) 



Eq. (|Al|l assumes that the density of states is a constant; 
and this is correct as long as the energy (oj + -^—) is 
larger than the Coulomb gap, A. At energies less than 
A, the density of states at the chemical potential is not a 
constant, but instead has the form p(e) ~ \e\ d ~ 1 (e/e 2 ) d . 
The Coulomb gap is the energy at which the density 
of states reaches the value in the absence of Coulomb 
interaction; thus A ~ [e 2d ni mp /(e5)] 1 ^ d ~ 1 \ Using this 
density of states, we get, for A > -2— > u), an optical 
conductivity 

Re(aH) (x Lor 2 - d ~ ^^{w/lo). (A2) 

At a finite temperature, a(cu — 0) is finite (see above). If 
the temperature is low, (T,u) <§C the frequency de- 
pendent conductivity in Eq. (|A1|I has an extra Boltzmann 
factor, 

Re((7(w, T)) ~ cr(0, T) + 27re 2 n 1 2 mp (l - e~ ul ' T ) x 

(u;/eS 2 )^hy d ( W M. (A3) 

For high enough temperatures, T 3> (cj, ^^),or high 
enough frequencies, u > (T, ^-), Mott's result, Eq.QJ 
is obtained^ 

The above treatment assumes that there is no inelas- 
tic scattering (e.g., by phonons). At finite temperature, 
phonons (with characteristic frequency bj p h ~ 10 12 Hz) 
provide an additional relaxation mechanism. Electrons 
make transitions by emitting or absorbing phonons with 
energy of the order of ui re i ~ uj p he~ Xij ^ loc , so the main 
contribution to the optical conductivity from inelastic 
processes comes from frequencies of the order of w r ei- 
For such frequencies, we should use u> p h instead of w in 
Eq.||I3). 



Appendix B: EFFECTIVE ACTION 
CORRECTIONS FROM EXPANSION OF G Vn IN 



POWERS OF ip n 



We explain how corrections to the effective tunnel- 
ing action in Ea. (|34(l coming from the expansion of 
G Vil (xi, Xj; r, r') in powers of ipn are small compared to 
the bare value when T>5. We expand Eq. (|2(J|) , 



where 



G i (0) (x i ,x; ; r,r')=T^ 



A,n 



V , A(x-)^( X c -^»(r-r') 
iv n - £iA 



= ^G i ( ° ) (^r')VA(xO^(x i ), 
I 

up to second order in ipn. The resulting correction to the 
effective action of Eq. (|35|l is 



t- 



AmR 2 



E 



xG^(r 1 ,rOG^ A2 (r',r)(n u+a (r,r'))^ 2 1 (r 1 )+ 

ia;AiA2A3 1 - 

xG. ( 2 (n, r 2 )G iAl (r 2 , r')G^ a ^ 3 (r',r)(n M+a (r, r')) x 

XVil( T l) ' Vil( T 2), 

(Bl) 

and Ili i_|_ a (r, r') is as defined in Ea. H36|) . Next we sim- 
plify Ea. l)Bl|l by integrating over r, r' . Using the fre- 
quency representation and completing the integration 
over t, t' we have 

SS&fai] = E / <n i ,i + aK l ))^ 2 1 (r 1 )x 

ia;Ai A2 ;n,m Tl 

1 1 



|<| 2 T 2 



(*% - €iAi) 2 [i{Vn + W m ) - 6+a,A 2 ] 

E / Gg ) 2 (n,T 2 ) ¥ » il (r 1 ). Vu (r 2 )x 

A 2 A 3 ;n,m ,/T1 ' r2 
e-fen^-rt) |p? AiA2 | 2 (n Li+a (c m )) 



2m 2 i? 2 , , 

ia;AiA2A3;n,m 1 



{iv n ~ £iAi) 2 [i{v n + W m ) - £i+a,A 3 ] ' 

(B2) 

It is convenient to perform the Matsubara sum over the 
fermionic frequencies. We have 



_. p-Wn(T2-Tl) 



1 



{iV n - ii\) 2 [i{v n + LO m ) - £i+a,A'] 



ia> TO + ^iA — £i+a,A' 



(B3) 



The first term in Eq. (|B2() vanishes when we use Ea. (IB3l) 
with ti =T2- Hence 



\t\ 2 T 



2m 2 R 2 E J b°AiAa 
ia;AiA 2 A 3 ;m T1T2 



G i V(T 1 ,T 2 Vii(Ti) ■Vii(r 2 )(n I ,i +a (a; m )): 



9 Gffi(7S,r 1 )-Gi"i 1 .(7i,r 1 )e 



(0) 



iw m (T 2 -ri) 
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Ea. (IB4|) contains two terms: one where the Green func- 
tions are on the same grain and the other where they are 
on different grains. The term with the Green functions 
on the same grain can be simplified by summing over 
£i+a,A 3 - The result of the summation in is(ep)sgn u> m . 
Since (Ili^+a^m)) is an even function of uj m , summing 
over io m makes the first term disappear. Thus so far, 



\t?T 



E / <$>(™) 



2m 2 B? 

>(0) 



x IVn(ri) • ¥>ii(r 2 )G^ a>A3 (r 2) n ) 

Q gio) m (T 2 — -ri) 



(n M+a (w m )) 



(B5) 



Now we integrate Ea. (|B5|) by parts with respect to the 
variable £i Al using J d£,i\i an d the identity 

d/d^x, = 2m^ A ( 1 /Pi,AA 1 )(9/5p i] A 1 A), 



WTv{e F ) 



•»v . ^ J d ^ J T < iT1,T2)x 



X<pn (n) • <pn (r 2 )G[+ ) a! A3 (r 2 , ti ) x 



x(n M+a (w m ))- 



iw m (r 2 -ri) 



iw m + £iAi — £i+a,A 3 

Now we complete the integration over £i Al to get 
2ir\t\ 2 Tiy(e F ) 



(B6) 



1R 2 



E / G ^M^ 



ia;A2A 3 ;m ,/T i' r2 

■r(O) 



X (pil (n ) • Vil (T-2 )G 1 Va, A 3 ( T 2 » T l ) X 

x (n iji+a (a; m ))sgn (uj m ) sin(w m (r 2 - n)). 

(B7) 



The right hand side of Eq. i|B7|l vanishes because the in- 
tegrand is odd with respect to interchange of t% and t 2 . 
This proves that the correction to the effective action 
arising from <pn fluctuations in G Vil may be ignored. 



Appendix C: OPTICAL CONDUCTIVITY OF 
THE AES MODEL 



We calculate the paramagnetic and diamagnetic terms 
in the AES conductivity ctaes of the granular array. As 
a corollary, we also find that the resonance width T is 
proportional to the diamagnetic part of the AES con- 
ductivity. From Eq.ljniJ through Eq.Jf^J it follows that 



1 



Re(a A Es(^T)) = -lm(K^ s (u;,T) + K^{u,T)). 

(CI) 

In the expression for Kj£g s (u),T) we need to calculate 
the cosine correlator, II(r) = (ITij+e^ (r)). We discuss 
the case of weak intergrain tunneling j< 1 first. 



(a) For weak intergrain tunneling, II (t) may be eval- 
uated perturbatively in increasing powers of g : H(t) = 
II( )(t) + Tv- 1 '(t) + ■ ■ ■ , where the prefixes denote the 
power of g. The leading term n(°)(r) can be shown to 
be2i 



nW(r) 



1 

Z 2 



E 



-/9Bc(«?+«i)-(l-?l-?2)2B T 



91,92 = — oo 



(C2) 



We similarly expand the diamagnetic response function 
Kf° s (m n ) in powers of g, K^ s = K d ^ s W + K d ™f ] + 
• • • , where the prefixes in brackets denote the power of g. 
To obtain the leading order in g behavior we use Eq. I|C2|I 
in Eq. l|62|l and take the Fourier transform: 



K 



dia,(l) 
AES ' 



ifl n ) 



( e V«)^EEi^i e 



-<?( 



z 2 



: ( e -2B ^(i-9i-92) _ J dnS(n-2E c (l-q 1 -q 2 ))x 

1 1 



(C3) 



Next we perform the Matsubara sum over m followed by 
an analytical continuation i£l n — ► u>. The result is 

Im(<^ ( W , T)) = g(e 2 /a)± £ e~^+^ x 



dfl(l 



Z 2 

n- 



9192 



coth ■ 



Cl- 



xS(n-2E c {l- qi -q 2 )). 



coth — 
2 

(C4) 



The DC (u> = 0) behavior in Ea. (|C4|) is dominated by 
single-charge excitations, (<?i, g 2 ) = (1,0), (0, 1), whereas 
the a.c. behavior at T = is dominated by "even" exci- 
tations {qi,q 2 ) = (0,0), (1, 1) : 

Im(^ s (1) (0,T))«2c 3 (e 2 / a )e-^ /T , 
Im(A^ a s (1 V,0)) « cog(e 2 /a)(l - 2E c /\u\)x 

X Q(\u\ - 2E C ). (C5) 

In the next order in <?, the cosine correlator can be 
shown 24 to be 



n«(r) 



2^T 2 



E 2 sin 2 (itTt) 



E c t > 1. 



(C6) 



It follows that in the second order in g, the imaginary 
part of the spectral function K 



dia,(2) 
AES 



IS 



Im(A^ s (2) (^,T)) = - 



47r.g 2 T 2 e 2 
3aE 2 



■[T 2 + (uj/2tt) 2 }, u^.E c . 

(C7) 



The power law behavior of the second order (in g) dia- 
magnetic response does not mean that the conductivity 
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will follow a power law. This is because we also have a 
paramagnetic contribution, and one can show<24 that the 
leading order paramagnetic response is second order in g 
and is equal and opposite to ^aeS > 



A AES V) ~ A 



para, (2) , 
AES I 



-A' 



dia,(2) 
AES 



(r). 



(C8) 



Thus power law contributions cancel out in the conduc- 
tivity and we are left with 



(C9) 



CAES 1,^7 1 ) > 

U) 

where Im(A^g (1) (w, T)) is given by Eq.£3|). 



Appendix D: BEHAVIOR OF THE RESONANCE 
WIDTH 



Now we discuss the frequency and temperature depen- 
dence of the polarization resonance width T. Note that 
T(u),T) is proportional to af^ s (uj,T). In the absence of 
a canceling paramagnetic contribution, T(u>,T), unlike 
the conductivity, does show an algebraic behavior at low 
frequencies, 



T(w,T) w Abzg- 



47re 2 
3RE? 



[T 2 + {uj/2-k)%uj <^E C . (Dl) 



Consider finally the case where the dimensionless in- 
tcrgrain tunneling is large, j > 1. Except at very low 
temperatures (explained below), both ctaes and T are 
more or less determined by the diamagnetic contribution. 
Evaluating the cosine correlator, 



n(r) » 1 HgE c r), gE c r » 1, 

ngz 



and substituting in the expression for diamagnetic re- 
sponse, we have 



ctaes(^,T) ps g{e 2 /a) 



i-J-mf *f - 

7T(7z \max(u;,Tj 



(D2) 



At exponentially small (in g) temperatures, such that the 
two terms in the square brackets in Eq. I|D2() become com- 
parable, perturbation theory in 1/g breaks down. Below 
such small temperatures, the behavior of oaes (and Y) is 
the same as for the g <C 1 case, except that the charging 
energy E c in g <C 1 results should now be replaced with 
an effective charging energyS&Si E*(g) that is exponen- 
tially small (in g) compared with E c . 
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